Near chromosome-level and highly repetitive genome assembly of the snake pipefish Entelurus aequoreus (Syngnathiformes: Syngnathidae).

The snake pipefish, Entelurus aequoreus (Linnaeus, 1758), is a northern Atlantic fish inhabiting open seagrass environments that recently expanded its distribution range. Here, we present a highly contiguous, near chromosome-scale genome of E. aequoreus. The final assembly spans 1.6 Gbp in 7,391 scaffolds, with a scaffold N50 of 62.3 Mbp and L50 of 12. The 28 largest scaffolds (>21 Mbp) span 89.7% of the assembly length. A BUSCO completeness score of 94.1% and a mapping rate above 98% suggest a high assembly completeness. Repetitive elements cover 74.93% of the genome, one of the highest proportions identified in vertebrates. Our demographic modeling identified a peak in population size during the last interglacial period, suggesting the species might benefit from warmer water conditions. Our updated snake pipefish assembly is essential for future analyses of the morphological and molecular changes unique to the Syngnathidae.

Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid?

Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?

The de novo genome of the Black-necked Snakefly (Venustoraphidia nigricollis Albarda, 1891): A resource to study the evolution of living fossils.

Snakeflies (Raphidioptera) are the smallest order of holometabolous insects that have kept their distinct and name-giving appearance since the Mesozoic, probably since the Jurassic, and possibly even since their emergence in the Carboniferous, more than 300 million years ago. Despite their interesting nature and numerous publications on their morphology, taxonomy, systematics, and biogeography, snakeflies have never received much attention from the general public, and only a few studies were devoted to their molecular biology. Due to this lack of molecular data, it is therefore unknown, if the conserved morphological nature of these living fossils translates to conserved genomic structures. Here, we present the first genome of the species and of the entire order of Raphidioptera. The final genome assembly has a total length of 669 Mbp and reached a high continuity with an N50 of 5.07 Mbp. Further quality controls also indicate a high completeness and no meaningful contamination. The newly generated data was used in a large-scaled phylogenetic analysis of snakeflies using shared orthologous sequences. Quartet score and gene concordance analyses revealed high amounts of conflicting signals within this group that might speak for substantial incomplete lineage sorting and introgression after their presumed re-radiation after the asteroid impact 66 million years ago. Overall, this reference genome will be a door-opening dataset for many future research applications, and we demonstrated its utility in a phylogenetic analysis that provides new insights into the evolution of this group of living fossils.

Genomic analysis reveals limited hybridization among three giraffe species in Kenya.

BACKGROUND: In the speciation continuum, the strength of reproductive isolation varies, and species boundaries are blurred by gene flow. Interbreeding among giraffe (Giraffa spp.) in captivity is known, and anecdotal reports of natural hybrids exist. In Kenya, Nubian (G. camelopardalis camelopardalis), reticulated (G. reticulata), and Masai giraffe sensu stricto (G. tippelskirchi tippelskirchi) are parapatric, and thus, the country might be a melting pot for these taxa. We analyzed 128 genomes of wild giraffe, 113 newly sequenced, representing these three taxa. RESULTS: We found varying levels of Nubian ancestry in 13 reticulated giraffe sampled across the Laikipia Plateau most likely reflecting historical gene flow between these two lineages. Although comparatively weaker signs of ancestral gene flow and potential mitochondrial introgression from reticulated into Masai giraffe were also detected, estimated admixture levels between these two lineages are minimal. Importantly, contemporary gene flow between East African giraffe lineages was not statistically significant. Effective population sizes have declined since the Late Pleistocene, more severely for Nubian and reticulated giraffe. CONCLUSIONS: Despite historically hybridizing, these three giraffe lineages have maintained their overall genomic integrity suggesting effective reproductive isolation, consistent with the previous classification of giraffe into four species.

The genome of the pygmy right whale illuminates the evolution of rorquals.

BACKGROUND: Baleen whales are a clade of gigantic and highly specialized marine mammals. Their genomes have been used to investigate their complex evolutionary history and to decipher the molecular mechanisms that allowed them to reach these dimensions. However, many unanswered questions remain, especially about the early radiation of rorquals and how cancer resistance interplays with their huge number of cells. The pygmy right whale is the smallest and most elusive among the baleen whales. It reaches only a fraction of the body length compared to its relatives and it is the only living member of an otherwise extinct family. This placement makes the pygmy right whale genome an interesting target to update the complex phylogenetic past of baleen whales, because it splits up an otherwise long branch that leads to the radiation of rorquals. Apart from that, genomic data of this species might help to investigate cancer resistance in large whales, since these mechanisms are not as important for the pygmy right whale as in other giant rorquals and right whales. RESULTS: Here, we present a first de novo genome of the species and test its potential in phylogenomics and cancer research. To do so, we constructed a multi-species coalescent tree from fragments of a whole-genome alignment and quantified the amount of introgression in the early evolution of rorquals. Furthermore, a genome-wide comparison of selection rates between large and small-bodied baleen whales revealed a small set of conserved candidate genes with potential connections to cancer resistance. CONCLUSIONS: Our results suggest that the evolution of rorquals is best described as a hard polytomy with a rapid radiation and high levels of introgression. The lack of shared positive selected genes between different large-bodied whale species supports a previously proposed convergent evolution of gigantism and hence cancer resistance in baleen whales.

Genomics reveals broad hybridization in deeply divergent Palearctic grass and water snakes (Natrix spp.).

Understanding speciation is one of the cornerstones of biological diversity research. Currently, speciation is often understood as a continuous process of divergence that continues until genetic or other incompatibilities minimize or prevent interbreeding. The Palearctic snake genus Natrix is an ideal group to study speciation, as it comprises taxa representing distinct stages of the speciation process, ranging from widely interbreeding parapatric taxa through parapatric species with very limited gene flow in narrow hybrid zones to widely sympatric species. To understand the evolution of reproductive isolation through time, we have sequenced the genomes of all five species within this genus and two additional subspecies. We used both long-read and short-read methods to sequence and de-novo-assemble two high-quality genomes (Natrix h. helvetica, Natrix n. natrix) to their 1.7 Gb length with a contig N50 of 4.6 Mbp and 1.5 Mbp, respectively, and used these as references to assemble the remaining short-read-based genomes. Our phylogenomic analyses yielded a well-supported dated phylogeny and evidence for a surprisingly complex history of interspecific gene flow, including between widely sympatric species. Furthermore, evidence for gene flow was also found for currently allopatric species pairs. Genetic exchange among these well-defined, distinct, and several million-year-old reptile species emphasizes that speciation and maintenance of species distinctness can occur despite continued genetic exchange.

Seascape Genomics and Phylogeography of the Sailfish (Istiophorus platypterus).

Permeable phylogeographic barriers characterize the vast open ocean, boosting gene flow and counteracting population differentiation and speciation of widely distributed and migratory species. However, many widely distributed species consists of distinct populations throughout their distribution, evidencing that our understanding of how the marine environment triggers population and species divergence are insufficient. The sailfish is a circumtropical and highly migratory billfish that inhabits warm and productive areas. Despite its ecological and socioeconomic importance as a predator and fishery resource, the species is threatened by overfishing, requiring innovative approaches to improve their management and conservation status. Thus, we presented a novel high-quality reference genome for the species and applied a seascape genomics approach to understand how marine environmental features may promote local adaptation and how it affects gene flow between populations. We delimit two populations between the Atlantic and Indo-Western Pacific oceans and detect outlier loci correlated with sea surface temperature, salinity, oxygen, and chlorophyll concentrations. However, the most significant explanatory factor that explains the differences between populations was isolation by distance. Despite recent population drops, the sailfish populations are not inbred. For billfishes in general, genome-wide heterozygosity was found to be relatively low compared to other marine fishes, evidencing the need to counteract overfishing effects. In addition, in a climate change scenario, management agencies must implement state-of-the-art sequencing methods, consider our findings in their management plans, and monitor genome-wide heterozygosity over time to improve sustainable fisheries and the long-term viability of its populations.

Range-wide whole-genome resequencing of the brown bear reveals drivers of intraspecies divergence.

Population-genomic studies can shed new light on the effect of past demographic processes on contemporary population structure. We reassessed phylogeographical patterns of a classic model species of postglacial recolonisation, the brown bear (Ursus arctos), using a range-wide resequencing dataset of 128 nuclear genomes. In sharp contrast to the erratic geographical distribution of mtDNA and Y-chromosomal haplotypes, autosomal and X-chromosomal multi-locus datasets indicate that brown bear population structure is largely explained by recent population connectivity. Multispecies coalescent based analyses reveal cases where mtDNA haplotype sharing between distant populations, such as between Iberian and southern Scandinavian bears, likely results from incomplete lineage sorting, not from ancestral population structure (i.e., postglacial recolonisation). However, we also argue, using forward-in-time simulations, that gene flow and recombination can rapidly erase genomic evidence of former population structure (such as an ancestral population in Beringia), while this signal is retained by Y-chromosomal and mtDNA, albeit likely distorted. We further suggest that if gene flow is male-mediated, the information loss proceeds faster in autosomes than in X chromosomes. Our findings emphasise that contemporary autosomal genetic structure may reflect recent population dynamics rather than postglacial recolonisation routes, which could contribute to mtDNA and Y-chromosomal discordances.

A chromosome-scale reference genome assembly of the great sand eel, Hyperoplus lanceolatus.

Despite increasing sequencing efforts, numerous fish families still lack a reference genome, which complicates genetic research. One such understudied family is the sand lances (Ammodytidae, literally: "sand burrower"), a globally distributed clade of over 30 fish species that tend to avoid tidal currents by burrowing into the sand. Here, we present the first annotated chromosome-level genome assembly of the great sand eel (Hyperoplus lanceolatus). The genome assembly was generated using Oxford Nanopore Technologies long sequencing reads and Illumina short reads for polishing. The final assembly has a total length of 808.5 Mbp, of which 97.1% were anchored into 24 chromosome-scale scaffolds using proximity-ligation scaffolding. It is highly contiguous with a scaffold and contig N50 of 33.7 and 31.3 Mbp, respectively, and has a BUSCO completeness score of 96.9%. The presented genome assembly is a valuable resource for future studies of sand lances, as this family is of great ecological and commercial importance and may also contribute to studies aiming to resolve the suprafamiliar taxonomy of bony fishes.

A Chromosome-Scale Genome Assembly of the Okapi (Okapia johnstoni).

The okapi (Okapia johnstoni), or forest giraffe, is the only species in its genus and the only extant sister group of the giraffe within the family Giraffidae. The species is one of the remaining large vertebrates surrounded by mystery because of its elusive behavior as well as the armed conflicts in the region where it occurs, making it difficult to study. Deforestation puts the okapi under constant anthropogenic pressure, and it is currently listed as "Endangered" on the IUCN Red List. Here, we present the first annotated de novo okapi genome assembly based on PacBio continuous long reads, polished with short reads, and anchored into chromosome-scale scaffolds using Hi-C proximity ligation sequencing. The final assembly (TBG_Okapi_asm_v1) has a length of 2.39 Gbp, of which 98% are represented by 28 scaffolds > 3.9 Mbp. The contig N50 of 61 Mbp and scaffold N50 of 102 Mbp, together with a BUSCO score of 94.7%, and 23 412 annotated genes, underline the high quality of the assembly. This chromosome-scale genome assembly is a valuable resource for future conservation of the species and comparative genomic studies among the giraffids and other ruminants.

Genomic Impact of Whaling in North Atlantic Fin Whales.

It is generally recognized that large-scale whaling in the 19th and 20th century led to a substantial reduction of the size of many cetacean populations, particularly those of the baleen whales (Mysticeti). The impact of these operations on genomic diversity of one of the most hunted whales, the fin whale (Balaenoptera physalus), has remained largely unaddressed because of the paucity of adequate samples and the limitation of applicable techniques. Here, we have examined the effect of whaling on the North Atlantic fin whale based on genomes of 51 individuals from Icelandic waters, representing three temporally separated intervals, 1989, 2009 and 2018 and provide a reference genome for the species. Demographic models suggest a noticeable drop of the effective population size of the North Atlantic fin whale around a century ago. The present results suggest that the genome-wide heterozygosity is not markedly reduced and has remained comparable with other baleen whale species. Similarly, there are no signs of apparent inbreeding, as measured by the proportion of long runs of homozygosity, or of a distinctively increased mutational load, as measured by the amount of putative deleterious mutations. Compared with other baleen whales, the North Atlantic fin whale appears to be less affected by anthropogenic influences than other whales such as the North Atlantic right whale, consistent with the presence of long runs of homozygosity and higher levels of mutational load in an otherwise more heterozygous genome. Thus, genome-wide assessments of other species and populations are essential for future, more specific, conservation efforts.

Speciation and population divergence in a mutualistic seed dispersing bird.

Bird-mediated seed dispersal is crucial for the regeneration and viability of ecosystems, often resulting in complex mutualistic species networks. Yet, how this mutualism drives the evolution of seed dispersing birds is still poorly understood. In the present study we combine whole genome re-sequencing analyses and morphometric data to assess the evolutionary processes that shaped the diversification of the Eurasian nutcracker (Nucifraga), a seed disperser known for its mutualism with pines (Pinus). Our results show that the divergence and phylogeographic patterns of nutcrackers resemble those of other non-mutualistic passerine birds and suggest that their early diversification was shaped by similar biogeographic and climatic processes. The limited variation in foraging traits indicates that local adaptation to pines likely played a minor role. Our study shows that close mutualistic relationships between bird and plant species might not necessarily act as a primary driver of evolution and diversification in resource-specialized birds.

Gene losses in the common vampire bat illuminate molecular adaptations to blood feeding.

Vampire bats are the only mammals that feed exclusively on blood. To uncover genomic changes associated with this dietary adaptation, we generated a haplotype-resolved genome of the common vampire bat and screened 27 bat species for genes that were specifically lost in the vampire bat lineage. We found previously unknown gene losses that relate to reduced insulin secretion (FFAR1 and SLC30A8), limited glycogen stores (PPP1R3E), and a unique gastric physiology (CTSE). Other gene losses likely reflect the biased nutrient composition (ERN2 and CTRL) and distinct pathogen diversity of blood (RNASE7) and predict the complete lack of cone-based vision in these strictly nocturnal bats (PDE6H and PDE6C). Notably, REP15 loss likely helped vampire bats adapt to high dietary iron levels by enhancing iron excretion, and the loss of CYP39A1 could have contributed to their exceptional cognitive abilities. These findings enhance our understanding of vampire bat biology and the genomic underpinnings of adaptations to blood feeding.

Conservation Genomics of Two Threatened Subspecies of Northern Giraffe: The West African and the Kordofan Giraffe.

Three of the four species of giraffe are threatened, particularly the northern giraffe (Giraffa camelopardalis), which collectively have the smallest known wild population estimates. Among the three subspecies of the northern giraffe, the West African giraffe (Giraffa camelopardalis peralta) had declined to 49 individuals by 1996 and only recovered due to conservation efforts undertaken in the past 25 years, while the Kordofan giraffe (Giraffa camelopardalis antiquorum) remains at <2300 individuals distributed in small, isolated populations over a large geographical range in Central Africa. These combined factors could lead to genetically depauperated populations. We analyzed 119 mitochondrial sequences and 26 whole genomes of northern giraffe individuals to investigate their population structure and assess the recent demographic history and current genomic diversity of West African and Kordofan giraffe. Phylogenetic and population structure analyses separate the three subspecies of northern giraffe and suggest genetic differentiation between populations from eastern and western areas of the Kordofan giraffe's range. Both West African and Kordofan giraffe show a gradual decline in effective population size over the last 10 ka and have moderate genome-wide heterozygosity compared to other giraffe species. Recent inbreeding levels are higher in the West African giraffe and in Kordofan giraffe from Garamba National Park, Democratic Republic of Congo. Although numbers for both West African and some populations of Kordofan giraffe have increased in recent years, the threat of habitat loss, climate change impacts, and illegal hunting persists. Thus, future conservation actions should consider close genetic monitoring of populations to detect and, where practical, counteract negative trends that might develop.

Whole-genome analysis of giraffe supports four distinct species.

Species is the fundamental taxonomic unit in biology and its delimitation has implications for conservation. In giraffe (Giraffa spp.), multiple taxonomic classifications have been proposed since the early 1900s.1 However, one species with nine subspecies has been generally accepted,2 likely due to limited in-depth assessments, subspecies hybridizing in captivity,3,4 and anecdotal reports of hybrids in the wild.5 Giraffe taxonomy received new attention after population genetic studies using traditional genetic markers suggested at least four species.6,7 This view has been met with controversy,8 setting the stage for debate.9,10 Genomics is significantly enhancing our understanding of biodiversity and speciation relative to traditional genetic approaches and thus has important implications for species delineation and conservation.11 We present a high-quality de novo genome assembly of the critically endangered Kordofan giraffe (G. camelopardalis antiquorum)12 and a comprehensive whole-genome analysis of 50 giraffe representing all traditionally recognized subspecies. Population structure and phylogenomic analyses support four separately evolving giraffe lineages, which diverged 230-370 ka ago. These lineages underwent distinct demographic histories and show different levels of heterozygosity and inbreeding. Our results strengthen previous findings of limited gene flow and admixture among putative giraffe species6,7,9 and establish a genomic foundation for recognizing four species and seven subspecies, the latter of which should be considered as evolutionary significant units. Achieving a consensus over the number of species and subspecies in giraffe is essential for adequately assessing their threat level and will improve conservation efforts for these iconic taxa.

Two high-quality de novo genomes from single ethanol-preserved specimens of tiny metazoans (Collembola).

BACKGROUND: Genome sequencing of all known eukaryotes on Earth promises unprecedented advances in biological sciences and in biodiversity-related applied fields such as environmental management and natural product research. Advances in long-read DNA sequencing make it feasible to generate high-quality genomes for many non-genetic model species. However, long-read sequencing today relies on sizable quantities of high-quality, high molecular weight DNA, which is mostly obtained from fresh tissues. This is a challenge for biodiversity genomics of most metazoan species, which are tiny and need to be preserved immediately after collection. Here we present de novo genomes of 2 species of submillimeter Collembola. For each, we prepared the sequencing library from high molecular weight DNA extracted from a single specimen and using a novel ultra-low input protocol from Pacific Biosciences. This protocol requires a DNA input of only 5 ng, permitted by a whole-genome amplification step. RESULTS: The 2 assembled genomes have N50 values >5.5 and 8.5 Mb, respectively, and both contain ∼96% of BUSCO genes. Thus, they are highly contiguous and complete. The genomes are supported by an integrative taxonomy approach including placement in a genome-based phylogeny of Collembola and designation of a neotype for 1 of the species. Higher heterozygosity values are recorded in the more mobile species. Both species are devoid of the biosynthetic pathway for ß-lactam antibiotics known in several Collembola, confirming the tight correlation of antibiotic synthesis with the species way of life. CONCLUSIONS: It is now possible to generate high-quality genomes from single specimens of minute, field-preserved metazoans, exceeding the minimum contig N50 (1 Mb) required by the Earth BioGenome Project.

A new genus and species of mud snake from Myanmar (Reptilia, Squamata, Homalopsidae).

Based on two male and two female individuals, we describe a new genus and species of mud snake, Myanophis thanlyinensis gen. nov., sp. nov., from the vicinity of the campus of East Yangon University, Yangon, Thanlyin, Myanmar. This species differs from every other homalopsid species by the following combination of characters: (1) dorsal scales smooth, row formula 21-21-19 or 21-21-17; (2) tail short, ratio tail length/SVL 0.185-0.204 in males, 0.160-0.167 in females; (3) nasal scales separated; (4) 125-126 ventral scales in males, 120-122 in females; (5) 38-39 subcaudal scales in males, 32-34 in females; and (6) hemipenis bilobed. Its matrilineal genealogy (based on analyses of 16S and cytochrome b sequences), associates Myanophis thanlyinensis gen. nov., sp. nov. most closely with species of the genera Myrrophis and Gyiophis. The new taxon differs from the species of Myrrophis and Gyiophis by having a bilobed hemipenis (vs. unilobed). Myanophis thanlyinensis gen. nov., sp. nov. differs further from the species of Myrrophis by having 125-126 ventral scales in males and 120-122 in females (vs. 137-162 and 137-164, respectively), and 38-39 subcaudal scales in males and 32-34 in females (vs. 39-55 and 37-52, respectively). Myanophis thanlyinensis gen. nov., sp. nov. differs further from the species of Gyiophis by lacking dark blotches along flank (vs. present), and by having 21 dorsal scales rows at midbody (vs. 25). We provide an identification key to the homalopsid species known to occur in Myanmar. As a novelty to the classic holotype description and characterization, the individual has been genome sequenced by Illumina short-read technology and its genome has been assembled into a draft nuclear genome and a complete, annotated mitochondrial genome. This innovative approach comprehensively and permanently characterizes the genomic variation of the holotype.

SambaR: An R package for fast, easy and reproducible population-genetic analyses of biallelic SNP data sets.

SNP data sets can be used to infer a wealth of information about natural populations, including information about their structure, genetic diversity, and the presence of loci under selection. However, SNP data analysis can be a time-consuming and challenging process, not in the least because at present many different software packages are needed to execute and depict the wide variety of mainstream population-genetic analyses. Here, we present SambaR, an integrative and user-friendly R package which automates and simplifies quality control and population-genetic analyses of biallelic SNP data sets. SambaR allows users to perform mainstream population-genetic analyses and to generate a wide variety of ready to publish graphs with a minimum number of commands (less than 10). These wrapper commands call functions of existing packages (including adegenet, ape, LEA, poppr, pcadapt and StAMPP) as well as new tools uniquely implemented in SambaR. We tested SambaR on online available SNP data sets and found that SambaR can process data sets of over 100,000 SNPs and hundreds of individuals within hours, given sufficient computing power. Newly developed tools implemented in SambaR facilitate optimization of filter settings, objective interpretation of ordination analyses, enhance comparability of diversity estimates from reduced representation library SNP data sets, and generate reduced SNP panels and structure-like plots with Bayesian population assignment probabilities. SambaR facilitates rapid population genetic analyses on biallelic SNP data sets by removing three major time sinks: file handling, software learning, and data plotting. In addition, SambaR provides a convenient platform for SNP data storage and management, as well as several new utilities, including guidance in setting appropriate data filters. The SambaR source script, manual and example data set are distributed through GitHub: https://github.com/mennodejong1986/SambaR.